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ABSTRACT 

The Pipe Nebula, a large nearby molecular cloud lacks obvious signposts of star formation 
in all but one of more than 130 dust extinction cores that have been identified within it. In 
order to quantitatively determine the current level of star formation activity in the Pipe Nebula, 
we analyzed 13 square degrees of sensitive mid-infrared maps of the entire cloud, obtained with 
the Multiband Imaging Photometer for Spitzer (MIPS) at wavelengths of 24 /im and 70 /im to 
search for candidate Young Stellar Objects (YSOs) in the high-extinction regions. We argue that 
our search is complete for class I and typical class II YSOs with luminosities of Lboi ^ 0.2 Lq 
and greater. We find only 18 candidate YSOs in the high-extinction regions of the entire Pipe 
cloud. Twelve of these sources are previously known members of a small cluster associated with 
Barnard 59, the largest and most massive dense core in the cloud. With only six candidate 
class I and class II YSOs detected towards extinction cores outside of this cluster, our findings 
emphatically confirm the notion of an extremely low level of star formation activity in the Pipe 
Nebula. The resulting star formation efficiency for the entire cloud mass is only ^ 0.06%. 

Subject headings: stars: formation — stars: pre-main sequence — infrared: stars 



1. Introduction 

The most critical but least understood aspects 
of the star formation process concern its initial 
conditions and earliest stages. In particular, little 
is known about the initial phases of protostellar 
development within dense molecular cloud cores 
and even less is known about the origin of the 
dense cores themselves. This is largely due to an 
apparent absence of nearby {d < 2.5 kpc) massive 



(> 10^ Mq) molecular clouds in a sufficiently early 
evolutionary state so that little or no star forma- 
tion activity is present within them. This is then 
taken as evidence that the onset of star formation 
in molecular clouds is extremely rapid and that 
clouds in the earliest evolutionary states must be 
exceedingly rare. Clearly the identification and 
detailed study of such clouds provides invaluable 
information about the earliest phases of star for- 
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mation. 

With a mass of 10^ M p^ and a distance of 130 pc 
( Lombardi et al. 2006h . the Pipe Nebula is one of 
the massive clouds nearest to the Sun and yet has 
received little attention and has only been stud- 
ied in detail within the last few years. With the 
exception of Barnard 59, the region is lacking the 
most obvious signposts of star formation, such as 
bright HII regions, reflection nebulae, luminous 
infrared sources and large populations of emission- 
line stars and may be a prime example of a molec- 
ular cloud in an early stage of evolutionary devel- 
opment. As such, it constitutes a critical labora- 
tory for s tar formation studies (for a review, see 
Alves et al.. .2008b.) . The first wide-area study of 
the region consisted of CO mapping carried out by 
Onishi et all (|l999 ') who identified 14 C^^O cores 
in a large low -resolution and unders ampled map. 
Subseauentlv. iLombardi et al. ( 20061 ) presented a 
high-resolution completely sampled dust extinc- 
tion map (based on 2MASS near-infrared data) to 
study general properties of the cloud. This extinc- 
tion map was also used to study its dense core pop- 
ulation: Alves et al. (2007) identified a population 
of roughly 150 extinction cores and found that the 
shape of the core mass function (CMF) is nearly 
identical to the stellar initial mass function (IMF) , 
but shifted to higher masses by a factor of three, 
suggesting that the CMF is a direct evolution- 
ary precursor to the st ellar IMF. Most recently. 



source (iReipurth et al. Il996h . and a CO outflow 



(|Onishi et al.lll999l ) . More recently, Spitzer obser- 



Rathborne et al 



(l2009l) have further refined the 
core identification by taking into account addi- 
tional measurements of C^^O strength and veloc- 
ity towards single cores as indicators of dense ma- 
terial and cloud membership. Their list contains 
134 dense cores in the Pipe Neb ula. To more reli- 
ably d etect the extinction peaks. [Rathborne et al 



()2009l) . lave used a background-subtracted version 
of the near-infrared extinction map, produced by 
filte ring out structures surpassing a given size (see 



also I Alves et al.llioOTh . Note that in the follow- 



ing, we we will generally refer to this background- 
subtracted map as the extinction map, except as 
otherwise noted. 

The only previously known region of active 
star formation in the Pi pe Nebu l a is the dark 
cloud Barnard 59 (B 59. barna^ Il927l ). Early 
studies found indications of ongoin g star forma- 
tion in B59: Ha emission line stars ( Herbid 20051 
and references therein), a millimeter continuum 



vations revealed a you ng cluster in side of B 59 



(|Brooke et all l2007l ). iRiaz et al.l (|2009l) stud 



led the cluster's most deeply embedded object, 
[BHB2007] 11, in more detail. In order to quan- 
tify the star formation activity for the whole cloud 
complex, we present a mid-infrared survey carried 
out with the Multiband Imaging Photometer for 
Spitzer (MIPS) as well as with the Infrared Array 
Camera (IRAC), both onboard the Spitzer Space 
Telescope. 

2. Observations and Data Reduction 

We present large-area mapping observations of 
the Pipe Nebula, carried out with Spitzer-MIPS 
(|Rieke et al.l 120041 ) . at wavelengths of 24 /im and 
70 /im. The observations were obtained between 
April 2004 and October 2006. The data consist 
of 44 Astronomical Observation Requests (AORs). 
They come from three different programs: the Cy- 
cle 1 GTO observations of B68 (Program Identifi- 
cation PID 53, 1 AOR), six mipsphot observations 
from Legacy observations (PID 139, c2d) and 35 
mipsscan as well as 3 mipsphot maps from the 
Cycle 2 GO Pipe Cloud Survey (PID 20119). We 
use Spitzcr-IRAC data that were mostly collected 
as part of program PID 20119, but also as pro- 
grams PID 139 and PID 30570. The observations 
were designed to cover the most highly extincted 
regions of the Pipe Nebula. In total, the 24 /im 
data cover an area of 13.4 square degrees, with the 
70 /im data covering a very similar area (to within 
2%). 

The data reduction was performed with the 
MOPEX software package, vlS.l.^l. The 24 /im 
data were self-calibrated by dividing each Basic 
Calibrated Data (BCD) image with a flatfield cre- 
ated by stacking all the BCDs in an AOR. All 
MIPS BCDs were from release S16. Self calibra- 
tion in this way removes large scale sky varia- 
tions and some artificats. The flat field image 
was created and applied using the fiatfield.pl tool 
provided as part of MOPEX. Briefiy, the BCD 
stack arc median filtered in order to detect sources, 
mask bright sources, and then median stack the 
original frames into a fiat, clipping out the bright 
sources. No overlap correction was applied to the 



^http: / /ssc. spitzer. caltcch.cdu/postbcd/mopcx. html 
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24 /xm data since flatfielding appeared to be suffi- 
cient. The AORs were individually processed into 
mosaics to preserve the two-epoch nature of most 
of the 24 /im data, allowing the rejection of as- 
teroids in the final source catalog. The mosaick- 
ing produced images with a final pixel size of 2.5" 
and included steps to mask faulty bits and reject 
outliers. No self-calibration or overlap corrections 
were applied to the 70 /im BCDs where we work 
with the median-filtered (fbcd) frames produced 
by the Spitzer Science Center pipeline process. 
The pixel size at 70 /im is 4". 

Source detection was performed on the individ- 
ual 24 /im mosaics, using MOPEX multi-frame 
PRF-fitting photometry. Given the Pipe Nebula's 
location close to the ecliptic and thus an increased 
probability to detect asteroids, we required source 
detections in both of the two observation epochs, 
if two observing epochs were available. The con- 
catenated list of source detections contains nearly 
80000 sources above a signal-to-noise ratio of four. 
In contrast, the 70 /im map mostly has only single- 
epoch coverage. Source detection was carried out 
on a single, large map produced from the median- 
filtered BCD files. In this map, a total of 5308 
sources were detected above a signal-to-noise ratio 
of four. A cross-correlation of the two lists with 
a search radius of 4" yields a list of 546 sources 
detected in both MIPS bands. 

For the identified candidate YSOs, we obtained 
Spitzer-IRAC mid-infrared aperture photometry 
from corresponding post-basic calibrated data 
(PBCD) sets for use in the fitting of their spectral 
energy distributions (SEDs), allowing a photo- 
metric accuracy of ^ 10%. Post-BCD mosaics 
from the Spitzer Science Center (pipeline S14) 
were used without modification. Aperture radii 
used were either two or three pixels (the former 
only for sources 22 and 24 to exclude neighboring 
sources), with a pixel size of 1.2" and the cor- 
responding aperture corrections. In a few cases, 
this required the use of a tooH to rectify saturated 
sources. 

While we processed and analyzed the crowded 
young cluster region Barnard 59 in the same way 
as all other data, we m ake use of the IRAC and 
MIPS photometry from iBrooke et"all (|2007l ) for 
two reasons: first, the crowding is especially prob- 



lematic in the MIPS data that we use for source 
detection, causing us to miss a few sources. How- 
ever, no comparable crowding occurs in the en- 
tire remaining areas in the Pipe Nebula. Sec- 
ond, the most clearly saturated sources are in 
B 59, requiring more effort to extrac t reliable pho- 
tomet ry. As was pointed out by iBrooke et al 



(2001), the three brightest 24 /im sources in B 59 
([BHB2007] 1, 7, and 11) are saturated while at 
70 /im, only [BHB2007] 11 is saturated. Outside 
of B 59, among the sources we are studying, only 
KK Oph is saturated at 24 /im. 

In addition to crowding of bright sources, our 
reliance on MIPS data for the identification of 
YSOs also suffers from the relatively low angular 
resolution: the detectability of multiple systems 
is reduced even when compared to 2MASS. For 
example, KK Oph is a known binary with a sep- 
aration of 1.6" (see below), but we cannot resolve 
the two components. 

3. Analysis 

3.1. Selection of candidate YSOs 

We expect to find any candidate protostars in 
the relatively high-extinction regions of the cloud. 
These regions correspond to only a small per- 
centage of t he full coverage. In line with ear- 
lier analysis ( Rathborne et al.1l2009l in particular), 
we use the following threshold to define the high- 
extinction regions: in an extinction map in which 
extended structure h as been removed b y means of 
wavelet subtraction ( Alves et al. 2007( ) we use a 



minimum extinction of Ay > 1.2 mag. An area of 
0.98 square degrees is above this threshold out of 



which the identified dense cores (Rathborne et al 



^ http: / / spider . ipac. caltech. edu / staff / j arrett /irac/ tools / 



2009t ) comprise an area of 0.75 square degrees. In 
the following, we will constrain the star forma- 
tion activity in these high-extinction regions. It 
will become clear that we can constrain the sizes 
of the populations of both cla ss I protost ars and 
class II YSOs (as defined by iLadal 119871 ) in the 
high-extinction regions. 

Given the location of the Pipe Nebula, pro- 
jected onto the vicinity of the Galactic center, the 
background source density in the region is enor- 
mous and selection criteria for candidate YSOs 
are of particular importance. Compared to source 
detection at near- infrared wavelengths, the lower 
source density in the mid-infrared at 24 /im and 
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70 /Ltm allows us to work with a sensitivity-limited 
(rather than confusion-limited) sample of objects. 
When looking for candidate protostars at these 
wavelengths, we first have to establish whether 
the mid-infrared mosaics are sensitive enough for 
this purpose. A good first test is to look at the 
lowest-luminosity protostellar objects known, the 
Very Low Luminosity Objects (VeLLOs). These 
objects were only discovered recently, with the un- 
precedented sensitivity of Spitzer. The four VeL- 



2004 



Bourke et all 120061 : 



in detail (lYoun 


E et al. 


Dunham et al. 


20061 



Lee et al.ll2009l . , L1014-IRS, L1521F-IRS, IRAM 
04191, and L328-IRS, respectively) have been de- 
tected with MIPS at both 24 /im and 70 /xm, and 
their mean bolometric luminosity is 0.25 Lq (the 
extreme values being 0.15 Lq, IRAM 04191, and 
0.36 Lq, L1521F-IRS). Note that the bolometric 
luminosities are significantly higher than the in- 
trinsic luminosities of the corresponding central 
objects. Detections at both 24 iim and 70 iim are 
also used as criteria bv lDu nham ct al" who 
list a total of 15 candidate VeLLOs as a result of 
surveying the Spitzer c2d data. The fact that the 
four above-mentioned VeLLOs are located nearby 
and at least as distant as the Pipe Nebula means 
that we can use them to define flux cut-off values 
at both 24 /zm and 70 /zm to include these cur- 
rently known VeLLOs. These flux limits ensure 
that we can detect even the lowest-mass protostel- 
lar objects. 

As our first selection criterion we thus require 
candidate YSOs to have flux densities of at least 
10 mJy at 24 /zm and 100 mJy at 70 fim. It turns 
out that both cut-off values are comfortably above 
the detection limits in both bands. In the list of 
sources detected in both MIPS bands, the typical 
photometric uncertainties determined by MOPEX 
are 0.09 mJy at the 24 /im flux cutoff and 5.65 mJy 
at the 70 fim flux cutoff. As a second step, we 
are using a direct application o f the original def- 



(Il987h . al- 



inition of YSO SED classes of iLada 
though implemented in longer-wavelength (MIPS) 
bands. Thus, we require a positive spectral index, 
a = dlogXFx/dlogX, in the two MIPS bands, in or- 
der to find the earhest (class I) protostellar stages. 
In this case, a positive spectral index in XF\ corre- 
sponds to a flux ratio of F,y(70 /im)/F^(24/im) > 
70/24. A similar criterion can be used to select 
candidate class II sources, dominated by their disk 



emission. They are expected to have flux ratios of 
F,y(70/im)/F^(24/im) > 1/3, corresponding to a 
spectral index of a > — 2. Sources above our flux 
cut-off values but below this minimum spectral 
index include candidate class III sources. Since 
these presumably become difficult to distinguish 
from stellar photospheres of field stars, we do not 
discuss these in further detail, focusing only on 
candidate class I and II sources as possible prod- 
ucts of recent star formation. Only these sources, 
as defined by our MlPS-based criteria, will be re- 
ferred to in the following as candidate YSOs. 

Using these criteria, we identify in Fig.[l]a total 
of only 25 candidate YSOs in the high-extinction 
regions of the entire complex, 12 of which are lo- 
cated in B 59. These sources were matched with 
the 2MASS point source catalog, using a search 
radius of 3". The median difference in position is 
1.5". We found one additional candidate YSO in 
a region without coverage at 70 /im. This source 
is listed as candidate 26. Lacking 70 /im cov- 
erage, source 26 was identified by its Ks — [24] 
color of 10.459±0.082 mag. There are no further 
sources with comparably extreme Ks — [24] colors 
that are not already covered by the above crite- 
ria. Given the superior sensitivity of the above 
method, which does not depend on near-infrared 
detections, we do not further discuss Ks — [24] col- 
ors. The resulting 26 candidate YSOs and their 
photometry are listed in Tables [T] and [21 sorted 
by descending 70/im flux. The photometric errors 
consist of the fitting errors reported by MOPEX, 
with the absolute error^ added in quadrature (for 
the saturated sources, a larger error is assumed). 
In order to allow us to compare the photometry, 
we have added these ab solute uncerta i nties also 
to the errors reported bv lBrooke erall (|2007t ) for 
the non-saturated sources, since the errors in their 
source table are the formal fitting errors. For all 
candidate YSOs we confirmed their detection in 
two different epochs of MIPS24 observations, en- 
suring that they were not caused by asteroids. 
Source 26 was confirmed by the presence of a mid- 
infrared IRAC counterpart. 

For B59, wh en applying o u r crit eria to the 



12 sources from iBrooke et al.l (|2007n that have 



MIPS24 and MIPS70 photometry, we find four 



^4% at 24 /^m and 7% at 70 fim, according to MIPS Data 
Handbook V3.3.1 



4 



candidate class I sources and eight candidate 
class II sources. Th is leaves eight sources from the 



Brooke et al.l (|2007l ) list of 20 that do not appear 



to qualify as candidate YSOs according to our 
MlPS-based criteria. While in some cases, these 
objects are genuinely excluded by our method, in 
other cases, no suitable MIPS photometry is avail- 
able. One of these sources, [BHB2007] 20 is not 
covered at 70 /zm so that it is not selected by our 
criteria. Three sources ([BHB2007] 2, 6, and 12) 
are too close to bright neighboring sources for any 
reliable photometry in either of the MIPS bands 
to be extracted (and n o such photometry is listed 
in Brooke et al.l l2007l ): therefore, these sources 
do not appear in our list. The remaining four 
sources are the outermost objects which are not 
de tected at 70 iiui; on ly upper limits were given 
bv lBrooke et"al] (|2007t l. These include [BHB2007] 
5 and 17, which Covey et al. (2009) determined 
to be background giants, using near-infrared spec- 
troscopy. The source [BIIB2007] 17 is very close 
to the chip edge in the 70 /Ltm coverage, and not 
detected in spite of its bright 24 counterpart. 

The sensitivity of our MIPS mosaics, as we have 
argued, is enough to even detect VeLLOs at the 
distance of the Pipe Nebula. Therefore, we es- 
timate that our survey is complete for candidate 
class I protostars of i > 0.2 Lq. We note that 
even for candidate class II sources, our selection 
method is very sensitive. Note that we used a con- 
servative limit of a24-7o > ^2 to find all possible 
candidate class II sources. This index corresponds 
to the ideal case of an optically thick, spatially 
flat disk. In reality, observed class II SEDs often 
have significantly fiatter slopes due to disk flaring 
and viscuous heating, result ing in a = —0.7 ± 0.3 
(jKenvon fc HartmanrJ 119871 for a spectral range 
of 2 - 60 iim). In Section [3^ when we study their 
SEDs, we will see from the final classification (Ta- 
bleHl) that our candidate class II sources have a 
mediar0 spectral index of aK~2i = —0.73. We 
can thus roughly estimate the completeness for 
such typical class II sources by comparing with 
an observed example, the fluxes of FM Tau, a K3 
class II source {aK-24 = —0.65) in Taurus at a 
similar distance as the Pipe. FM Tau has a to - 
tal luminosity of 0.5 L© (e.g., IStrom et al. 1988 ) 



and a 70 /zm flux density of 0.29 Jy (from the 
Taurus Spitzer Legacy Project). Scaling that to 
our flux density cut-off at 70 /im (0.1 Jy) and 
assuming the same distance for the sake of this 
comparison yields a completeness limit for class II 
SEDs of only Lboi ~ 0.2 Lq. The flux density ratio 
F^(70^m)/F^(24^m) in this case is 0.59. A sim- 
ilar check with two late-type class II objects, FN 
Tau (M5, aK-2i = -0.53) and FP Tau (M2.5, 
<^/f-24 — ^1-12), yields even lower limiting lu- 
minosities for our flux limits (i.e., 0.05 L© and 
0.12 Lq, respectively), corroborating our estimate 
of iboi ~ 0.2 Lq. The completeness level for typi- 
cal observed class II sources is comparable to the 
one for class I protostars discussed above. 

While all candidate sources by deflnition are in 
regions with Ay > 1.2 mag (in the extinction map 
without extended structure, se e above), not all 



of the m are projected onto what iRathborne et al 



(I2OO9I) define as dense cores, see Table [T] Five 
of our sources either lie in extinction cores that 
do not fulfill the m i nimum size requirement set 
bv IRathborne et ah ( 2009^ or were excluded due 
to their C^^O emission properties: One source is 
close to an extinction peak without detected C^'^^O 
emission (source 16), and, most interestingly, the 
extinction core onto which source 17 is projected, 
as well as its neigboring cores, have C^^O emission 
at a velocity that is significantly different from the 
cloud velocity {v = 21kms~'^ vs . a range for the 



*The values range from ctx— 24 = —0.54 to aK— 24, = —0.83, 
leaving aside source 24 with a markedly double-peaked SED 



(a 



K-2A 



-1.47). 



Pipe Nebula oi v — 2..8kms ^, IRathborne et al 
20091) . This candidate YSO therefore apparently 
does not belong to the Pipe Nebula complex. For 
the five cases where cores and t heir masses are 
not listed by IRathborne et all (|2009l ). Table [1] 
lists corresponding mass estimates, derived ac- 
cording to the procedure outlined in Section 13.61 
These estimates take into account all pixels of 
the background-subtracted extinction map with 
Ay > 1.2 mag if they have not been assigned to 
neighboring cores. Note that some of the changes 
in the refined core list are due to the intricacies of 
defining the core population, a process that is not 
unambiguous. 

3.2. Comparison to control regions 

While we have estimated the sensitivity of our 
selection method, we still have to determine the 
amount of contamination by background sources. 
As control region, we use those parts of our MIPS 
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mosaics that are not covering the high-extinction 
regions, i.e., regions with Ay < 1.2. Applying the 
selection criteria to these sources allows us to ap- 
proximately predict how many detections in each 
category could be due to background objects and 
to judge whether the number of identified candi- 
dates exceeds the number of expected background 
sources. As mentioned above, about 94% of our 
map is not covering high-extinction regions. Given 
the proximity of the Pipe Nebula to the Galac- 
tic center, there is a high density of background 
sources, presumably mainly in the Galactic bulge. 
While 56 sources were detected at both 24 /xm and 
70 fim in regions with extinctions above 1.2 mag 
(not applying any filtering) , 490 such sources were 
detected in the remainder of the map. Fig. [2] de- 
picts the sources in the control region in the same 
way as the high-extinction sources were plotted in 
Fig. [H Some of the sources in the control region 
technically fulfill our candidate YSO selection cri- 
teria, indicating that there is contamination in our 
list of candidate YSOs. Table[3]contains a compar- 
ison of the control region and the high-extinction 
regions in the Pipe Nebula in terms of candidate 
YSOs. Taking into account that regions with ex- 
tinctions above 1.2 mag cover only ~ 7% of the 
map (see above), we find that that in the three 
ranges of flux ratios corresponding to YSOs, we 
find more candidates than can be expected from 
background contamination alone. In particular, 
and most importantly, we detect nine candidate 
class I sources while the expected background con- 
tamination is 2.9 sources. This relative overden- 
sity of candidate YSOs in the high-extinction re- 
gions suggests that their occurence is indeed re- 
lated to the presence of dense material and sup- 
ports the notion that some of these objects really 
are YSOs. 

3.3. Comparison to other star-forming re- 
gions 

There are two obvious possibilities to compare 
the application of our candidate YSO identifica- 
tion method to other star-forming regions. First 
of all, given its location at about the same dis- 
tance, it is interesting to compare the Pipe Nebula 
to the Taurus star-forming region. While a full ap- 
plication of our method to the Taurus data is be- 
yond the scope of this paper, we can at least check 
whether known objects occupy the region of flux 



space that we use to identify candidate YSOs. In 
a preliminary version of their source list, the Tau- 
rus Spitzer legacy project^ lists 693 sources that 
are detected at both 24 /xm and 70 /zm, a pre- 
condition for our candidate YSO selection method. 
The Taurus map covers a total of about 40 square 
degrees, four times larger than our map of the 
Pipe Nebula. We correlate the Spitzer source list 
with the hst of YSO s in Taurus- Auriga given by 
Kenvon et al. Illooi). This table contains 383 low- 



mass YSOs of which 98 have counterparts in the 
Spitzer list with detections in both MIPS bands 
(within a search radius of 1", see above; note that 
the two lists cover different areas). The Kenyon et 
al. list does not distinguish the different evolution- 
ary stages, but as shown in Fig. [3l this selection 
of sources is well constrained by the regions in the 
70 iim vs. 24 /j,m plot that we use for the selection 
of candidate YSOs. To illustrate that prototypical 
class I sources fall into the corresponding selection 
region, we note that L1551 IRS5 (not covered by 
the Taurus Spitzer Legacy Project), already based 
on its IRAS data, fulfills the criteria, but lies out- 
side of our plot limits due to its brightness: its 
60/Ltm flux is 373 Jy, and its 25/xm flux is 106 Jy, 
i.e., the flux ratio is 3.5. This example once more 
underlines the sensitivity of the Spitzer MIPS ob- 
servations. 

For a second comparison, we consider the candi- 
date YSOs that have been identified by the "Cores 
to disks" (c2d) Spitzer Legacy Project. This com- 
parison is different, since it allows us to select al- 
ready identified candidate YSOs and look at their 
properties at 24 /im and 70 /im. The c2d "Ensem- 
ble Cores" list of candidate YSOs contains 106 ob- 
jects with both 24 /im and 70 /im fluxes out of 122 
candidate YSOs that were observed at both wave- 
lengthfl Reassuringly, when we plot the sources 
in the same way as we plotted the candidate YSOs 
in the Pipe, the sources roughly fall into the region 
of flux space where we expect candidate YSOs in 
the Pipe cloud (Fig. [3]). This is reassuring since 
the MIPS spectral index is not part of the selection 
criteria that were used by the c2d team. Note that 
the data have not been corrected for the different 



^ http: //data, spitzer . caltech.edu / popular / taurus2 / 

20081007.enlianced/ 
''Sixteen sources are listed as "undetected" at 70 fim, but 

are located in the vicinity of other, mostly bright, 70 /im 

sources, complicating their detection. 
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distances 120-250 pc) of the several different 
nearby star-forming regions that were part of the 
survey. 

These two examples demonstrate that using the 
24 ^m and 70 fim MIPS bands and their high 
sensitivity constitutes a powerful method to look 
for candidate YSOs in the Pipe cloud. The draw- 
back of relatively low angular resolution is in our 
case overcompensated by the enormous reduction 
of sources to survey (while keeping a high sensitiv- 
ity) compared to the use of IRAC data which in 
this direction of the Galaxy are affected by source 
crowding. 

3.4. A closer look at the candidate YSOs 

A cross-correlation of the 26 candidate YSOs 
with catalogs, using SIMBAD and Vizier, allows 
us to identify a first subset of sources that can 
be rejected as YSO candidates. Source 10 is a 
1.4 GHz radio continuum source in the NRAO 
VLA Sky Survey and hsted as a possible pre- 
planetarv nebula bv lKohoutek (2001). Sources 13, 
14, and 15 are all listed as hi gh mass-loss asymp - 
totic giant branch stars by lOiha et al.l ( 20071 ). 



After a comparison with typical Spitzer galaxy 
SEDs, in particular looking at the influence of 
PAH features on the IRAC photometry (e.g., 



Rowan-Robinson et al. 20051 ). three additional 



sources, namely 19, 20, and 22 have been re- 
jected as being likely galaxies. We note that 
their IRAC colors are [5.8/^m]-[8.0/im]>2 and 
[3.6 jUm] - [4.5iUm]<0.9, suppor ting this conclusion 
(e.g., iRobitaille et all l2008h . Source 25 looks 
photospheric except for its 70 /im flux which is 
the weakest in this selection. In fact, the 70 ^m 
detection is surrounded by weak extended emis- 
sion leading to an overestimate of the flux in the 
PRF fit. The source has been rejected as not be- 
ing a true excess source; it is clearly not a class 
I or class II YSO. The use of catalogs also allows 
us to confirm the selection of a YSO: source 5 is 
the Herbig AeBe star KK Oph which is very likely 
associated t o B 59; it ap p ears to have a T Ta. u 
companion (|Herbigl [iooj ICarmona et"all l2007^ . 
Additionally, Covey et al. (2009) have analyzed 
near-infrared spectroscopy of our sources 11 and 
16, confirming only the latter as a YSO (their 
source 2, a binary system unresolved in our data), 
while leaving open the identific a,tion of the former 
(source 1 in their designation). iRiaz et "aL I (|2009l) 



confirm the protostellar nature of source 1. Thus, 
after this first cross-correlation with catalogs and 
previously obtained data as well as comparison 
to galaxy SEDs, 18 candidate YSOs remain, only 
six of them outside of B 59. We note that the two 
iBrooke et al. (|2007l ) sources rejected as giant stars 
by Covey et al. (2009), [BHB2007] 5 and 17, are 
not detected at 70 /im and thus did not appear in 
our list. 

In order to better characterize the candidate 
YSOs, we have determined their infrared SEDs, 
based on 2MASS, IRAC and MIPS photometry. 
For comparison, we have also determined SEDs 
for some sources in the control regions that would 
qualify as candidate YSOs. Subsequently, to test 
whether the SEDs can be explained as protostel- 
lar, we have attempted to fit th e SEDs of can- 



didate YSOs with models from Robitaille et al 



(|2007n . These pre-computed radiative-transfer 
models cover a large region of parameter space. 
Since most of the sources are relatively bright 
2MASS sources and thus do not look like edge- 
on disk configurations, we excluded such models 
in the grid from the fitting process, realizing that 
all SEDs can still be well fitted and explained. For 
the fits, photometric errors larger than the nomi- 
nal errors were assumed: 0.1 mag for 2MASS and 
10% for all Spitzer photometry except for data 
corrected for saturation, where we assume 20%. 
The SEDs and fit results are shown in Figs. [H 
[51 and [HI Table 2] lists the range of total ex- 
tinction and the range of bolometric luminosities 
for all candidate YSOs. The parameter ranges 
are simply giving the extreme values for all fits 
within a relatively strict goodness-of-fit critcrior0, 
~ Xhest < '^data- The fit rcsults show that the 
candidate YSOs can be reasonably described as 
a mostly low-mass YSO population. The lowest- 
luminosity candidate class I protostar, source 26, 
might be a VeLLO in its own right. Also, the low- 
est luminosities deduced from the model fits are 
of the same order as our estimated completeness 
limit. 

We use the full SEDs and the model fits to 
reexamine and refine our initial classification of 
the candidate YSOs which was based solely on 



^For source 9, only a single fit was found with this criterion; 
we therefore relax the criterion for only this source to — 
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the MIPS spectral index. In particular, Figs. |4] 
and [5] show that two candidate YSOs, sources 9 
and 24, initially identified as candidate class I 
sources, exhibit complex double-peaked SEDs that 
cannot be adequately described by a single spec- 
tral index. The full SEDs show, however, that 
these two sources are more aptly described as 
class II objects. Additionally, an examination of 
the SEDs of sources initially selected as candidate 
class II objects shows that several of these have 
SEDs that are flatter than expected for prototyp- 
ical class II objects with shapes similar to those 
of the fiat-sp e ctrum protostellar class defined by 
Greene et~aU (|l994[ ). Since their SEDs are well- 
behaved, we use the wide K-24 spectral index, 
as described in a footnote to Tabled to refine our 
classification to include flat-spectrum protostars. 
The refined YSO classifications are listed in Ta- 
bleU 

In an attempt to characterize the 48 control 
sources that technically fulfill our criteria for be- 
ing candidate class I protostars, we used SIMBAD 
to identify counterparts. Such counterparts are 
found for only nine sources within a search radius 
of 4". Eight of these appear to be planetary neb- 
ulae and one is an OH/IR star. These 48 sources 
have a mean extinction of 0.1 mag in the extinction 
map where extended structure has been removed, 
and a mean extinction of 4.0 mag in the full ex- 
tinction map. A sample of control SEDs is shown 
in Fig. [SI While many sources remain unidenti- 
fied and the multi-wavelength coverage is not the 
same for all of these objects, most of the sampled 
sources appear to be galaxies based on their SED 
shapes, following the same arguments also used in 
rejecting candidate YSOs above (mainly the influ- 
ence of PAH features on the IRAC photometry). 

3.5. Location of the candidate YSOs rela- 
tive to extinction peaks 

Very young candidate YSOs would most likely 
be found at or in close proximity of the extinc- 
tion peaks. Fig. [7] shows the identified candi- 
date YSOs overlaid onto the corresponding sec- 
tions of the extinction map. Given that the can- 
didate YSOs were selected for being inside areas 
of the extinction map that are above a minimum 
threshold, it is not surprising to see them close 
to the extinction peaks, but few are centered on 
peaks, the best case being source 26. Since one 



would expect the youngest sources (class I) to be 
located at or very near to the extinction peaks, 
the observed offsets indicate that the candidate 
YSOs are likely predominantly in later evolution- 
ary stages (i.e., class II). This is compatible with 
our flux-based characterization which found that 
most of the candidate YSOs are class II rather 
than class I (see above). It is interesting to note 
that while five of the six candidate YSOs outside of 
B 59 are within 2.3° of the cluster, a single candi- 
date (source 17) is 5.9° away, at the "opposite end" 
of the cloud. We note again that the correspond- 
ing core as well as the surroundi ng cores were re- 
jected bv lRathborne et aL I (I2OO9I ) as being part of 
the Pipe complex since their C^^O velocities are 
significantly offset from the system velocity. While 
source 17 nevertheless is a candidate YSO, it prob- 
ably does not belong to the Pipe Nebula complex. 
Fig. [8] shows the six candidate YSOs outside of 



B 59 o n the extinction map from iLombardi et al 
( 20061 ). It is indeed interesting to note that most 



of them are located in the "stem" of the Pipe, none 
having been detected in the "bowl" , even though 
there are regions of considerable extinction in that 
area. This suggests that the area around B 59 at 
the tip of the stem is the most evolved region, fol- 
lowed by the stem, leaving the entire "bowl" in a 
very early evolutionary stage. Interestingly, this is 
in line with the conclusions of Alves et al.( ( 2008a ) 
who presented an optical polarization study of the 
Pipe Nebula. From a strong polarization gradient 
across the Pipe Nebula, with the lowest amounts 
of polarization found towards B 59, they propose a 
very similar picture. While they also suggest that 
the entire "bowl" region is in a very early stage, 
they conclude that the "stem" is about to form 
stars while B 59 has actively ongoing star forma- 
tion. 



3.6. Star formation efficiency in the Pipe 
Nebula 

In total, we have found 21 candidate YSOs in 
the Pipe Nebula. Fifteer[f| of these are in B 59 
and only six in the entire remaining cloud com- 
plex. Star formation efhciencics (SFEs) are usu- 



*In addition to the 12 MlPS-selected candidate YSOs, we 
consider the three above-mentioned iBrooke et al. ] ll2007f) 
sources in B 59 for which no photometry could be obtained 
due to nearby bright neighbors, yielding an upper limit of 
15 YSOs in B59. 
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ally given as SFE= M^/{M^ + Mdoud), i-e., they 
are compared to the total mass of t he cloud (in this 
case, M ~ 10^ Mf^. lOnishi et all [1999). Assum- 
ing stellar masses of ~ 0.3 Mq each, this trans- 
lates into an extremely low star formation effi- 
ciency when compared to the mass of the entire 
cloud, it is only ~ 0.06%. Typical val ues are a few 



it is much higher when compared to the dense 
material , as derived from dust continuum obser- 
vations ( Evans et al. 20091 ). 



percen t; the median of 2% given by iMvers et al 



( 19861 ) is 33 times higher than our value for the 
Pipe Nebula. 

Furthermore, we can estimate the star forma- 
tion efficiency for the densest material, consider- 
ing that the high-extinction regions represent the 
dense component of the cloud. We can estimate 
the corresponding cloud mass from the summed 
pixels in the extinction map, taking into account 
the high-extinction regions with Ay > 1.2 mag 
(in the background-corrected map) that were sur- 
veyed for candidate YSOs: 



M = /^mH/3^pixSAy 



(1) 



Here, /Lt is the mean molecular weight (1.37), 
mn is the mass of the hydrogen atom, f3 is the 
conversion from optical extinction to column den- 
sity, and T,Av is the summed optical extinction. 
Apix is the area of a pixel in cm~^ at the distance 
of the source, for extinction pixels with a side 
length of spix ["]: ^pix = (71/180/3600)2^;^^^ 
where d is the distance in cm. While there are 
14066 pixels in the extinction map above a thresh- 
old of Av > 1.2 mag (each of size 30" x 30"), 
only 11063 of these are covered by both the 24 /Ltm 
and the 70 /im maps. The summed Ay in these 
pixels is 29751 mag. Using a conversion factor o f 
f3 = 2x10^1 cm-^mag-i fe.g- lBohlin et al.lll978l ). 
the resulting mass is 230 Mq. Assuming an upper 
limit of 21 candidate YSOs, each approximately of 
^ 0.3 Mq, again yields an extremely low star for- 
mation efficiency of 2.7% when compared to the 
densest material. This is considerably lower than 
typical values for clusters which are o n the order 
of 10-40% (e.g., iLada fc Ladalll99lh . We note 
that at most five of the identified 134 extinction 
cores (or 3.7%) have candidate YSOs, not count- 
ing three additional candidate YSOs in extinction 
cores projected onto hig h-extinction regions no t 
part of cores as defined bv lRathborne et al ( 20091 ). 
For the combined data of the Spitzer c2d project, 
the average star formation efficiency is 4.8% while 



4. Summary 

We have used very sensitive mid-infrared map- 
ping of the Pipe Nebula to obtain a sensitive cen- 
sus of its extremely low star formation activity. 
Using the 24 /^m and 70 /im MIPS bands to iden- 
tify candidate YSOs is a powerful tool, allowing us 
to search the cloud for protostars down to objects 
with very low luminosities. For both class I and 
typical class II YSOs the search is estimated to be 
complete down to luminosities of ^ 0.2 Lq. The 
star formation efficiency in the Pipe Nebula is only 
~ 0.06% when compared to the total cloud mass 
and 2.7% when compared to the mass of the high- 
extinction regions. This is more than an order of 
magnitude lower than in typical star-forming re- 
gions. Candidate YSOs are found in at most five 
of the 134 identified extinction cores, leaving the 
overwhelming majority of cores apparently star- 
less. Outside of Barnard 59, where a previously 
known cluster of YSOs resides, we find only six 
new candidate YSOs in the entire remaining Pipe 
Nebula, one of which has already been confirmed 
as a YSO by near-infrared spectroscopy (Covey et 
al. 2009, in prep.). One of these candidate YSOs 
very likely is not a member of the Pipe Nebula 
complex. Using our MlPS-based selection crite- 
ria, we find an upper limit of 15 candidate YSOs in 
B59, all of them previously known. Interestingly, 
the Pipe candidate YSOs outside of B 59 are not 
evenly distributed over the remaining cloud, but 
they are within 2.3° of B59, leaving large regions 
of the cloud, particularly the bowl region, without 
any candidate YSOs. These results clearly demon- 
strate that apart from B 59, the Pipe Nebula is a 
region with extremely low star formation activity. 
We conclude that it appears to be an ideal region 
to study the initial conditions and earliest stages 
of star formation. 

We would like to thank Thomas Robitaille for 
help with the SED fits. This work is based on ob- 
servations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology under a 
contract with NASA. Support for this work was 
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University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/ California Institute 
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Table 1 

Candidate young stellar objects in the Pipe Nebula 



no. 




RA, Dec 


(J2000)'' 


core 


core mass 






MIPS70 


err. 


Apes'* 

L J 


MIPS24 


err. 


i="„(70 Aim)/ 


ID" 












[ Mn 1 


[niag] 


[mag] 


[Jy] 






[Jy] 


F„ (24 fira) 




1 


17: 


11:23.09 - 


-27:24:32.80 


9 


19.4 


20.2 


24.9 


42 


,100 


8.400 


0.5 


4.715 





,940 


8.9 


[BHB2007] 11 


2 


17: 


11:17.45 - 


-27:25:08.48 


9 


19.4 


11.1 


15.0 


10 


,100 


0.709 


3.6 


9.200 


1 


,840 


1.1 


[BHB2007] 7 


3 


17: 


11:22.12 - 


-27:26:02.30 


9 


19.4 


8.8 


9.9 


9 


,190 


0.644 


0.7 


1.710 





,074 


5.4 


[BHB2007] 10 


4 


17: 


11:03.95 - 


-27:22:55.39 


7 


4.0 


4.2 


7.9 


8, 


,490 


0.596 


0.6 


6.320 


1 


,260 


1.3 


[BHB2007] l' 


5 


17: 


10:08.20 - 


-27:15:20.20 




<0.1 


2.4 


3.0 


2 


,901 


0.203 


1.5 


6.100 


1 


,250 


0.5 


KK Oph 


6 


17: 


11:21.43 - 


-27:27:42.91 


9 


19.4 


12.4 


15.7 


2, 


,070 


0.145 


1.4 


1.540 





,064 


1.3 


[BHB2007] 9 


7 


17: 


11:27.00 - 


-27:23:49.05 


9 


19.4 


10.5 


18.1 


1 


,780 


0.125 


0.6 


0.831 





,034 


2.1 


[BHB2007] 13 


8 


17: 


11:18.13 - 


-27:25:49.30 


9 


19.4 


7.5 


13.0 


1 


,440 


0.101 




0.412 





,017 


3.5 


[BHB2007] 8 


9 


17: 


11:30.32 - 


-27:26:29.57 


9 


19.4 


10.3 


14.4 


1 


,310 


0.092 


0.9 


0.150 





,006 


8.7 


[BHB2007] 16 


10 


17: 


35:02.22 - 


-25:23:49.54 




0.3 


2.2 


8.5 


0, 


,919 


0.065 


0.6 


2.160 





.086 


0.4 


R 


11 


17: 


14:51.05 - 


-27:22:41.02 


17 


1.6 


2.2 


5.6 


0, 


,817 


0.058 


0.8 


0.550 





,022 


1.5 


Src 1' 


12 


17: 


11:27.28 - 


-27:25:28.99 


9 


19.4 


12.4 


20.3 


0, 


,710 


0.050 


1.5 


0.455 





,018 


1.6 


[BHB2007] 14 


13 


17: 


:29:41.64 - 


-26:29:09.34 


49 


0.8 


2.1 


6.6 


0, 


,688 


0.049 


3.4 


0.906 





,036 


0.8 


R 


14 


17: 


:25:18.87 - 


-26:46:07.79 


40 


4.2 


1.7 


5.1 


0, 


,669 


0.047 


2.5 


1.237 





,049 


0.5 


R 


15 


17: 


:36:57.59 - 


-24:18:06.02 




0.3 


1.4 


4.9 


0, 


,469 


0.034 


0.5 


1.346 





,054 


0.3 


R 


16 


17: 


14:56.44 - 


-27:31:09.26 




0.5 


1.5 


4.1 


0, 


,323 


0.023 


0.2 


0.416 





.017 


0.8 


Src 2' 


17 


17: 


37:35.42 - 


-23:30:52.19 




1.5 


3.7 


7.9 


0, 


,270 


0.020 


1.3 


0.213 





.009 


1.3 




18 


17: 


11:29.31 - 


-27:25:36.30 


9 


19.4 


12.4 


15.6 





,246 


0.018 




0.104 





.004 


2.4 


[BHB2007] 15 


19 


17: 


:20:18.85 - 


-26:55:57.35 


29 


4.3 


1.4 


4.8 


0, 


,241 


0.017 


1.8 


0.014 





.001 


16.7 


R 


20 


17: 


:34:01.76 - 


-25:24:10.73 


92 


1.0 


1.7 


9.3 


0, 


,213 


0.015 


1.4 


0.013 





.001 


16.2 


R 


21 


17: 


11:14.45 - 


-27:26:54.40 


9 


19.4 


7.5 


11.8 


0, 


,210 


0.016 




0.155 





,006 


1.4 


[BHB2007] 4 


22 


17: 


23:35.15 - 


-23:42:24.05 


100 


2.4 


2.8 


4.1 


0, 


,184 


0.013 


0.8 


0.011 





,001 


16.8 


R 


23 


17: 


11:11.78 - 


-27:26:55.17 


9 


19.4 


9.0 


12.9 


0, 


,110 


0.009 


1.4 


0.095 





.004 


1.2 


[BHB2007] 3 


24 


17: 


:20:54.17 - 


-26:46:42.77 


34 


9.3 


2.3 


5.7 


0, 


,110 


0.009 


0.7 


0.014 





.001 


7.6 




25 


17: 


:29:28.24 - 


-25:56:25.08 


57 


2.8 


3.9 


8.6 


0, 


,104 


0.008 


2.6 


0.014 





.001 


7.6 


R 


26 


17: 


19:41.25 - 


-26:55:31.67 


29 


4.3 


8.1 


10.6 










0.163 


0, 


,007 







Note. — Sources sorted according to their 70 /J.m flux. Photometry for Barnard 59 sources is from lBrooke et al.l (120071) . 
^from MIPS 24 /im data except for sources 8, 18 and 21, which were added from IBrooke et"all I I2007I) , see text 

optical extinction from an extinction map with extended structure removed JAIves et al.ll2007D 
■^optical extinction from an extinction map including extended structure jLombardi et 311120061) 
"^distance between MIPS70 and MIPS24 positions in source detection, possibly misleading for bright sources 

— rejected candidate YSO (see text) 
^at MIPS wavelengths: blend with neighboring [BHB2007] 2; pair also known as LkHa 346 AB. 
^in Covey et al. 2009 (in prep.) 



Table 2 

Candidate young stellar objects in the Pipe Nebula: 2MASS counterparts 



no. 


zMAob name 


2MASS-J 


err. 


01\ T A CO U 

zMAob-ii 


err. 


2MASS-Ks 


err. 






[mag 


] 


[mag] 


[mag] 


[mag] 


[mag] 




[mag] 


1 


J17112317-2724315 


>18.78 




>17.80 




15 


.08 


0.14 


2 


J171 11726-2725081 


13.62 


0.04 


10.82 


0.04 


8, 


.77 


0.03 


3 
4 


J17110392-2722551 


10.46 


0.03 


8.99 


0.03 


7, 


.76 


0.02 


5 


J17100811-2715190 


9.07 


0.03 


7.23 


0.04 


5, 


.80 


0.03 


6 


J17112153-2727417 


12.74 


0.02 


10.57 


0.02 


8, 


.98 


0.03 


7 


J17112701-2723485 


11.88 


0.04 


10.14 


0.04 


9 


.08 


0.03 


8 


J17111827-2725491 


>18.4 




15.12 


0.08 


11 


.95 


0.03 


9 


J17113036-2726292 


11, 


,91 


0.03 


10.00 


0.03 


8, 


.89 


0.02 


10 


J17350225-2523503 


15. 


,55 


0.07 


14.34 


0.08 


13, 


.11 


0.06 


11 


J17145108-2722399 


10. 


,26 


0.02 


8.89 


0.02 


8, 


.16 


0.02 


12 


J17112729-2725283 


13. 


17 


0.03 


10.65 


0.03 


9 


.11 


0.02 


13 


J17294185-2629092 


13. 


,75 


0.04 


10.06 


0.02 


7, 


.63 


0.02 


14 


J17251886-2646099 


>16. 


,05 




13.23 


0.04 


9 


.72 


0.03 


15 


J17365757-2418066 


>12. 


,71 




10.65 


0.06 


8 


.40 


0.03 


16 


J17145648-2731080 


10. 


,42 


0.02 


9.56 


0.03 


9, 


.04 


0.02 


17 


J17373543-2330517 


14. 


,60 


0.05 


12.38 


0.04 


10, 


.92 


0.03 


18 


J17112942-2725367 


13. 


,28 


0.02 


11.74 


0.02 


10, 


.70 


0.02 


19 
20 


J17340187-2524070 


15. 


71 


0.08 


14.02 


0.07 


13, 


.31 


0.06 


21 


J17111445-2726543 


11. 


,62 


0.02 


10.42 


0.02 


9, 


.63 


0.02 


22 


J17233513-2342249 


16. 


10 


0.11 


15.06 


0.10 


14, 


.39 


0.10 


23 


J17111182-2726547 


14. 


,09 


0.03 


12.62 


0.03 


11 


.76 


0.03 


24 


J17205433-2646434 


12. 


,30 


0.02 


11.02 


0.02 


10, 


.46 


0.03 


25 


J17292824-2556252 


9, 


,01 


0.03 


7.93 


0.03 


7, 


.47 


0.03 


26 


J17194124-2655317 


>18, 


,35 




>15.86 




14 


.57 


0.08 


Note 


. — Sources 3 and 19 do not have 2MASS counterparts. 











Table 3 

Comparison of detections in the high-extinction regions and the control regions 





control region 

n 


scaled to high-ext. 

n 


area 


high-ext. 

n 


region 


"class I" 


48 




2.9 




9 


"class 11" 


64 




4.0 




16 


"class HI" 


87 




5.2 




9 



Note. — The designations in the first column correspond to the flux space 
regions shown in Figures fTI and [2] 
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Table 4 

Fit results for candidate young stellar objects in the Pipe Nebula 



Source Av(tot) Lhol ckk-24 class^ 

mag mag Lq Lq 

mill max min max 



1 


116, 


,7 


179, 


.9 


1, 


,4 


2, 


.6 


2, 


.71 


I 


2 


20, 


.7 


117, 


.7 


16, 


,5 


30, 


.1 


0, 


.58 


I 


3 


78, 


.2 


996, 


.9 


0, 


,7 


20, 


.1 






I 


4 


8, 


.1 


106, 


.3 


7, 


,0 


29, 


.6 


0, 


.03 


f 


5 


8, 


.7 


62, 


.9 


55, 


,9 


175, 


.4 


-0, 


,73 


II 


6 


14, 


.1 


24, 


.7 


2, 


,4 


11 


.5 


-0, 


,09 


f 


7 


8, 


.9 


12, 


.8 


0, 


7 


2, 


.7 


-0, 


,30 


f 


8 


40, 


.1 


48, 


.1 


1, 


,4 


12, 


.3 


0, 


.50 


I 


9 


9, 


.3 


13, 


.8 


0, 


5 


4, 


.5 


-1, 


,09 


II 


11 


4, 


.3 


34, 


.0 


0, 


,8 


6, 


.6 


-0, 


,83 


II 


12 


15, 


.3 


22, 


.4 


1, 


,0 


7, 


.6 


-0, 


,54 


II 


16 


1 


.9 


14, 


.8 


0, 


,6 


2, 


.5 


-0, 


,61 


II 


17 


15, 


.0 


16, 


.4 


0, 


,3 


0, 


.6 


-0, 


,17 


f 


18 


6, 


.6 


36, 


.7 


0, 


,1 


0, 


.8 


-0, 


,55 


II 


21 


4, 


.5 


21 


.9 


0, 


,2 


1 


.1 


-0, 


,79 


II 


23 


9, 


.5 


26, 


.4 


0, 


,1 


0, 


.2 


-0, 


,18 


f 


24 


1 


.1 


3, 


.2 


0, 


,1 


0, 


.2 


-1, 


,47 


II 


26 


30, 


.0 


31 


.7 


0, 


,3 


0, 


.3 


1 


.12 


I 



Note. — For source 9, the goodness-of-fit criterion has been re- 
laxed to allow for more than one fit (see text). 

^evolutionary class, determined from the full SEDs; aK-24 is 
used to differentiate flat-spectrum from class I and II sources 
dGreene et al.| 119941) : a > 0.3 ^ class 1; -0.3 < a < 0.3 
class f (fiat spectrum); and a < —0.3 — >■ class II 



Young, C. H., J0rgensen, J. K., Shirley, Y. L., 
et al. 2004, ApJS, 154, 396 



This 2-coIumn preprint was prepared with the AAS lAT^TjX 
macros v5.2. 
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Fig. 1. — Plot of the MIPS fluxes of sources in the high-exinction regions. Regions where we expect to find 
candidate Young Stellar Objects in various evolutionary stages are marked. The boundaries of these regions 
were adapted from the original definition of YSO classes at shorter wavelengths, as discussed in Section lSH 
Dashed error bars denote the four VeLLOs that have been studied in detail (see text). 



14 




Fig. 2. — Plot of the MIPS fluxes of sources in the control region of the maps, outside of the high-extinction 
regions {Ay < 1.2). The areas used for the identification of candidate Young Stellar Objects are marked, as 
in Fig. [1] even though we only use this plot to determine the background contamination. 
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Fig. 3. — Plot of the MIPS sources in other star forming regions: '+' symbols mark MIPS sources among 
the candidate YSOs from the ensemble list of cores produced by the Cores to disks Spitzer Legacy Project. 
The lines denote the same selection regions that we also used for the Pipe candidate YSOs. Note that the 
plot contains objects from star-forming regions at different distances, without a cor rection for this effect . 
Asterisks mark the MIPS fluxes of counterparts to YSOs in Taurus- Auriga, as listed bv lKenvon et al. ( 2008 ). 
see text. For comparison, the flux regions for our YSO selection are shown. Both the c2d sources as well as 
the Taurus sources mostly fall into the area that we use for candidate YSO identification. 
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Src 1 = [BHB2007] 11 



Src 2 = [BHB2007] 7 




Src 3 = [BHB2007] 10 



Src 4 = [BHB2007] 1 




1 10 100 1 10 100 

X C/im) X 

Src 6 = [BHB2007] 9 Src 7 = [BHB2007] 13 
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Fig. 4. — SEDs of candidate YSOs in B59. The black line shows the best fit with gray lines delineating 
alternative fits within the goodness-of-fit criterion used (see text). Note that for Source 9, this criterion had 
to be relaxed. 
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Src 8 = [BHB2007] 8 



Src 9 = [BHB2007] 16 




1 10 

A C/im) 

Src 21 = [BHB2007] 4 



1 10 
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Fig.H- cntd. 
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Src 5 = KK Oph Src ll^ 
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Fig. 5.— SEDs for candidate YSOs outside of B 59 
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Src 10 



Src 13 




Src 19 (note different scale) 



Src 20 




Fig. 6. — SEDs of rejected candidate YSOs and four sample control sources that fulfill the YSO selection 
criteria (designated only with their coordinates). 
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Fig. El - cntd. 
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-1.2 -1-3 -1.4 



Fig. 7. — Location of the identified candidate YSOs relative to extinction peaks. Extinction is shown in 
shades of grey, starting at Ay < 1.2 and advancing in steps of 0.5. Due to its high extinction, the B59 
region is instead shown in an Ay range from 2 to 20 in steps of 2. The coordinates are Galactic coordinates. 
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Fig. 8. — Location of the identified candidate YSOs on the full extinction map from lLombardi et al. ( 20061 ). 
Due to the high source density in B 59, only the location of the cluster is indicated. The approximate outline 
of the combined Spitzer-MIPS coverage (24 /im and 70 /im) is marked. For details, see Fig. [T] 
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